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LARGE DEVIATION RATE FUNCTIONS FOR THE PARTITION 
FUNCTION IN A LOG-GAMMA DISTRIBUTED RANDOM 

POTENTIAL 



O ■ NICOS GEORGIOU AND TIMO SEPPALAINEN 

^ [ Abstract. We study right tail large deviations of the logarithm of the partition function 

I for directed lattice paths in i.i.d. random potentials. The main purpose is the derivation of 

explicit formulas for the 1+1-dimensional exactly solvable case with log-gamma distributed 
■ random weights. Along the way we establish some regularity results for this rate function 

for general distributions in arbitrary dimensions. 



1. Introduction 



^' 

■ We study a version of the model called directed polymer in a random environment where 
a fluctuating path is coupled with a random environment. This model was introduced in 
the statistical physics literature in [16] and early mathematically rigorous work followed 

■ in [31 |T7] . We consider directed paths in the nonnegative orthant Z^J. of the d-dimensional 
I integer lattice. The paths are allowed nearest-neighbor steps oriented along the coordinate 

^ ■ axes. A random weight a;(u) is attached to each lattice point u G Z'^. Together the weights 

ly^ ■ form the environment io = {io{u) : u S Z+j- The space of environments is denoted by 17. P 

I is a probability measure on Q under which the weights {a;(u)} are i.i.d. random variables. 

For V, u € such that v < u (coordinatewise ordering) the set of admissible paths 
from V to u with |u — v|i = m is 

nv,u = {x. = {v = xq, xi, . . . , Xm = u} : VA;, x^ € Z'^ and 



Xfc+i - Xk e {ei : 1 < i < d}} 
where is the i-th standard basis vector of M'^. The point-to-point partition function is 
(1.2) Zv,,= eS^--(^'^). 

This is the normalization factor in the quenched polymer distribution 

m 

(1-3) Q.,ui^.) = Z-^^l[e-^^^^ 
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which is a probabihty distribution on the paths in the set IIv^u- When paths start at the 
origin (v = 0) we drop v from the notation; Zu = Zq^u and IIu = Hq.u- Note that the 
weight at the starting point xq was not included in the sum in the exponent in (ll.2p . This 
makes no difference for the results. Sometimes it is convenient to include this weight and 
then we write = e^^^^Z^^^ where the superscript □ reminds us that all weights in the 
rectangle are included. 

In the polymer model one typically studies fluctuations of the path and fluctuations of 
logZu- This paper considers only logZu- Specifically our main object of interest is the 
right tail large deviation rate function 

(1.4) Ju(r) = - lim n"MogP{log Z|„u| > nr} 

where u G W^, r G M and the floor of a vector is [nyj = ([nyij, [ny2ji ■ ■ ■ , ["'T-Ud])- This 
function J exists very generally for superadditivity reasons, and in Section [3] we establish 
some of its regularity properties. 

The focus of the paper is an exactly solvable case where d = 2 and — w(u) is log-gamma 
distributed. By "exactly solvable" we mean that special properties of the log-gamma case 
permit explicit computations, such as a formula for the limiting point-to-point free energy 

(1.5) p(y) = liiii n-^^^ogZinyi P — a.s. 

and fluctuation exponents [31]. In the same spirit, in this paper we compute explicit 
formulas for the rate function J and other related quantities in the context of the 1+1- 
dimensional log-gamma polymer. 

One can also consider point-to-line partition functions over all directed paths of a fixed 
length. For m € N the partition function is defined by 

(1.6) Z^-= Yl 

uGZ^J_:|u|i=m 

Due to the log in front, in the results we look at Zl^'^ behaves like the maximal Zu 
over |u|i = m. 

Some comments are in order. 

There are currently three known exactly solvable directed polymer models, all in 1+1 
dimensions: The two with a discrete aspect are (i) the log-gamma model introduced in jSlj , 
(ii) a model introduced in [27] where the random environment is a collection of Brownian 
motions. Some fluctuation exponents were derived for the second model in [32], and it has 
been further studied in [261 via a connection with the quantum Toda lattice. This Brownian 
model possesses structures similar to those in the log-gamma model, so we expect that the 
results of the present paper could be reproduced for the Brownian model. 

The third exactly solvable model is the continuum directed random polymer [T] that is 
expected to be a universal scaling limit for a large class of polymer models (see jlQ] for a 
recent review). 

Usually the directed lattice polymer model is placed in a space-time picture where the 
paths are oriented in the time direction. (See articles and lectures O [U [U [13] for recent 
results and reviews of the general case.) In two dimensions (1 time + 1 space dimension) 
the space-time picture is the same as our purely spatial picture, up to a 45° rotation of 
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the lattice and a change of lattice indices. The temporal aspect is not really present in our 
work. So we have not separated a time dimension but simply regard the paths as directed 
lattice paths. 

Another standard feature of directed polymers that we have omitted is the inverse tem- 
perature parameter (3 € (0, oo) that usually appears as a multiplicative constant in front 

of the weights: Zv,u = eHv u ^■'^P'f'^ Sj=i '^(^j)j'' ^ fixed weight distribution, (5 
modulates the strength of the coupling between the walk and the environment. It is known 
that in dimension 1+3 and higher there can be a phase transition. By contrast, in low 
dimensions (1+1 and 1+2) the model is in the so-called strong coupling regime for all 
< /? < oo [71 [21] . The (3 parameter plays no role in the present work and has a fixed value 
/3 = 1. This is the unique (3 value that turns the log- gamma model into an exactly solvable 
model. 

The techniques of the current paper are entirely probabilistic and rely on the stationary 
version of the log-gamma model. It can be expected that as a combinatorial approach 
to this model is fully developed [11], more complete results and alternative proofs for the 
present results can be found. 

Earlier literature. Precise large deviation rate functions for log Z in the case of directed 
polymers have not been derived in the past. The strongest concentration inequalities can 
be found in recent references [U [22| I33j . The normalization of the left tail varies with 
the distribution of the weights as demonstrated by [2] but the right tails have the same 
normalization n. [4J has some bounds on the left tail of log Z in Gaussian environments in 
dimensions 1 + 3 and higher and for small enough j3. Similar bounds were proved later in 
[24j for bounded environments using concentration inequalities for product measures. 

For the exactly solvable zero-temperature models (that is, last passage percolation mod- 
els) large deviation principles have been proved. For the longest increasing path among 
planar Poisson points, an LDP for the length resulted from a combination of articles 
[HI EOl [23l [30]. These results came before the advent of determinantal techniques. For 
the corner growth model with geometric and exponential weights [T8] derived an LDP in 
addition to the Tracy- Widom limit. An earlier right tail LDP appeared in [29| . 

Notation. We collect some notation and conventions here for easy reference. N is for 
positive integers, Z+ for nonnegative integer, M+ for nonnegative real numbers and is the 
set of all vectors with nonnegative real coordinates. Vector notation: elements of M'^ and 1/ 
are v = (ui, f2, . . . , Vd)- Coordinatewise ordering v < u means vi < ui,V2 < U2, ■ ■ ■ ,Vc[ < 
Ud- Particular vectors: 1 = (1, 1, . . . , 1) and = (0, 0, . . . , 0). [y\={[yi\,[y2\,---,[yd\) 
where [yj = max{n € Z : n < y} is the integer part of y € M. The norm on is 
|v|i = \vi\ H h \vd\. 

The convex dual of a function / : M — )• (—00,00] is f*{y) = sup^^^{xy — f{x)}, and 
f = f** iff f is convex and lower semicontinuous. We refer to [28] for basic convex analysis. 

The partition function Z does not include the weight of the initial point of the paths, 
while Z^ does. In 2 dimensions we write Zm,n = ■2^{m,n)- 

The usual gamma function is r(^) = x^~^e~^ dx for /Lt > 0. The digamma and 
trigamma functions are ^'o = F'/F and ^1 = ^'q. On (0, 00) ^'o is increasing and concave 
and ^1 decreasing, positive and convex, with — ^'o(0+) = ^'i(0+) = 00. 
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2. Large deviations for the log-gamma model 

2.1. The log-gamma model with i.i.d. weights. In this section we specialize to d = 2 
dimensions and the log-gamma distributed weights. Fix a positive real parameter /i. This 
parameter remains fixed through this entire section, and hence is omitted from most nota- 
tion. In the log-gamma case we prefer to switch to multiplicative variables. So the weight 
at point € is Yij = e^*-*'-^-* where the reciprocal has Gamma(//) distribution. 
Explicitly, 

/oo 
x^-^e-^' dx for s G M+. 
_ 

As above, we write Y for a generic random variable distributed as Yij. The digamma and 
trigamma functions give the mean and variance: E(logy) = — ^'o(^) and Var(logy) = 

The logarithmic moment generating function (l.m.g.f.) of a; = logl^ is 

riogr(^-0-iogr(/i), ^e{-oo,fi) 
\oo, ^ e [^,oo). 

The point-to-point partition function for directed paths from (0, 0) to (m, n) is 

m+n 

(2.3) Zm,n= Yl n^-r 

Note that we simplified notation by dropping the parentheses: Zm,n = Zf^^^^)- For {s,t) € 
the limiting free energy density exists by superadditivity: 

(2.4) p{s,t) = \\VD. log Z\ns\Ant\ IP-a.s. 
The limit is a finite constant. We begin by giving its exact value. 



(2.2) M^(0=logE(e«i°sn 



IS 



Theorem 2.1. For {s,t) G and /x G (0, oo), the limiting free energy density (12. 4p 
given by 

(2.5) p{s, t) = inf {-s^oip) - t^o(M - P)}. 

0<p</i 

The value p{s,t) was already derived in [31] but the proof was buried among estimates 
for fluctuation exponents. In Section U] we sketch an elementary approach that utilizes 
special features of the log-gamma model. For the other explicitly solvable 1-1-1 dimensional 
polymer with Brownian environment, [25] computed the limiting free energy with a very 
different large deviation approach. 

The next result is a large deviation principle (LDP) for log ^[nsj,[ntj under normalization 
n. The rate function is 

f sup {r(- mi {tMe{0-sM^_e(.-0)}, r>p{s,t) 

(2.6) /,,t(r) = <^ Ce[o,/.) sg^.m) 

[oo, r < p{s, t). 

On the boundary (s = or t = 0) the result reduces to i.i.d. large deviations so we only 
consider (s,t) in the interior of the quadrant. 
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Theorem 2.2. Let Y^^ ~ Gamma{ii) as in ()2.ip and {s,t) G (0,oo)^. Then the dis- 
tributions of n^^ log Z^ns\,lnt\ satisfy a LDP with normalization n and rate function Ls,t- 
Explicitly, these bounds hold for any open set G and any closed set F in M.- 

(2.7) Ih^^ n-ilogP{n-ilogZL„,j,L„tj G F} < - mf /s,t(r), 
and 

(2.8) lim n-ilogP{n-ilogZL„,j,LntJ G G} > - inf /,,t(r). 

n— )-oo reG 

On [p(s,t),oo) i/ie rate function Ls^t is finite, strictly increasing, continuous and convex. Ln 
particular, the unique zero of Is^t{f) is at r = p{s,t). The right tail rate defined in (|1.4p is 
given by 



(2.9) Js,tir) 



0, r G (-oo,p(s,t)] 

Is,ti'^)' ^ ^ b(s,t),oo). 



Remark 2.3. From a computational point of view, the solution to the variational problem 
in (|2.6p can be computed by 

Ls,t{r) = sup {fr{9) - inf fr{z)} = /.(^s) - fr{ei), 
o<e<ti o<z<e 

where 

fr{9) = r6 + tlogr{6) - slogVifi - 6), 

and for any r > p{s,t), < 9i < 62 < fJ. are the solutions to the equation ^/r(^) = 0. (See 
Fig{TJ) This again implies that the rate function is strictly positive as long as r > p{s,t). 

Remark 2.4. We do not address the precise large deviations in the left tail, that is, in the 
range r < p(s,t). We expect the correct normalization to be n^. (Personal communication 
from I. Ben-Ari.) Presently we do not have a technique for computing the rate function 
in that regime. We include the trivial part Is,ti'r') = 00 for r < p{s,t) in the theorem so 
that we can compute the limiting l.m.g.f. by a straightforward application of Varadhan's 
theorem. 

Define for ^ G M 
(2.10) KtiO= lim n"^logEe^'°s^L-J.LntJ. 

Corollary 2.5. Let ^ G M. Then the limit in (j2.10p exists and is given by 



(2.11) A,,i(0 = lUO = { 



inf UMeiO - sM^_ei-0}, < ^ < fi 
00, C ^ 
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Figure 1. Graphical representation of the solution to the variational prob- 
lem ()2.6p that gives the rate function J<s,j(r) = fr{&2) — /r(^i) • The curve 
fr{&) has the same general shape as long as r > p{s,t). 

Remark 2.6. Symmetry of Ag^t in (-5, t) is clear from (j2.10p but not immediately obvious 
in the < < ^ case of (|2.1ip . It turns out that if s < t the infimum is achieved at 
a unique 9q € + ^)/2,/i), and then for At^siO the same infimum is uniquely achieved 
at 6*1 = // + ^ — 6*0 G (C, (a* + 0/2]- In the case s = t a simple formula arises: Af^tiO = 

2i(iogr(^)-iogr(ii±«)). 

Remark 2.7. The first case of (|2.6p gives Is^t as the dual A*^, and the reader may wonder 
whether this is the logic of the proof of the LDP. It is not, for we have no direct way to 
compute Ag^f Instead, Theorem 12. 21 is first proved in an indirect manner via the stationary 
model described in the next subsection, and then A^^t is derived by Varadhan's theorem. 

Let us also record the result for the point-to-line case. It behaves like the point-to-point 
case along the diagonal. 

Corollary 2.8. Let Y^^ ~ Gamma{^) as in ()2.ip and s > 0. Then the distributions of 
logZj'^^^j satisfy a LDP with normalization n and rate function 1^/2, s/2- 

Remark 2.9. For e > and r = p{s,t) + e, one can show after some calculus that there 
exists a non zero constant C = Cs^tifJ-) so that 

I,,{r) = Ce'/' + o{ey'). 

This suggests that Var(log Z^^^j^^^^j ) is of order n^/^. Rigorous upper bounds on the mo- 
ments E| log Z^ns\,[nt\ — np{s,t)\P for 1 < p < 3/2 can be found in [3l], Theorem 2.4. 
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We computed the precise value of the constant C for the point-to-hne rate function, 

(2.12) hAr) = \ e"' + o{e'l% 

3^1^2^2)1 

where ^2 = ^o- 

2.2. The stationary log-gamma model. Next we consider the fog-gamma model in 
a stationary situation that is special to this choice of distribution. Working with the 
stationary case is the key to explicit computations, including all the previous results, and 
provides some explanation for the formulas that arose for I^^t and A<j^f in (j2.6p and (j2.1ip . 

The stationary model is created by appropriately altering the distributions of the weights 
on the boundaries of the quadrant Z^. We continue to use the parameter /i G (0, cxd) 
fixed at the beginning of this section, and we introduce a second parameter 9 G (0,^). 
Let the collection of independent weights {Uifl,Vo^j,Yij : i,j G N} have these marginal 
distributions: 

(2.13) Uj^Q ~ Gamma(^), Vq'^ ~ Gamma(/x — 9), and Y-~j^ ~ Gamma(/x). 

Define the partition function Zm'n by (j2.3p with the following weights: at the origin Yq^ = 1, 
on the X-axis Yi^ = Uifi, on the y-axis Yqj = Vqj, and in the bulk the weights {Yij : 
i,j G N} are i.i.d. Gamma(/i)~^ as before. Equivalently, we can decompose the stationary 
partition function Zm'n according to the exit point of the path from the boundary: 

m y k \ n , I \ 

(2.14) Z(?)„ = E n f^^-O ^g,l),(-,n) + E n ^0.^- ^Pl/),(™,n)- 

The symbols f/j^o and were at first introduced for the boundary weights to highlight 
the change of distribution. Next let us define for all (i,i) G \ {(0, 0)} 

(2.15) = and V,, = 

(0) 

Note that this property was already built into the boundaries because for example q = 
Uifl ■ ■ ■ Uifi. The key result that allows explicit calculations for this model is the following. 

Proposition 2.10. For each G Z\ \ {(0,0)} we have the following marginal distribu- 
tions: U^^ ~ Gamma{9) and Vfj^ ~ Gamma(fi — 9). For any fixed n G l^j^, the variables 
{Ui^n '■ i £ N} are i.i.d., and for any fixed m G Z+, the variables {Vmj '■ j G N} are i.i.d. 

This is a special case of Theorem 3.3 in ^1], where the independence of these weights 
along more general down-right lattice paths is established. Proposition 12.101 is the only 
result from |,31j that we use. It follows in an elementary fashion from the properties of the 
gamma distribution. 

As an immediate application we can write 

[ntj [nsJ 

(2.16) log zf2j,L.tj = n-' Vo, + n"^ ^ log C/,L„ij 

j=i i=i 
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as a sum of two sums of i.i.d. variables, and from this compute 

(2.17) E(log Z^)„) = mE(log U) + nE(log V) = -m^o(^) - ^^oIa* - 0) 
and obtain the law of large numbers: 

(2.18) n-MogZ[2j,LntJ -^P^'Hs,t) = -s*o(^) - t^oil^ - 0) P-a.s. 

Note that the two sums on the right of (j2.16p are not independent of each other. In fact, 
they are so strongly negatively correlated that the variance of their sum is of order ii?/^ 
|31j . Comparison of ()2.5p and (I2.18P reveals a variational principle at work: p{s,t) is the 
minimal free energy of a stationary system with bulk parameter fj,. 

Instead of the right tail large deviation rate function we give the asymptotic l.m.g.f. in 
the next result. Define 

(2.19) Ae(st)(C)= lim n-^logEe^^°^^L-J,LntJ. 

Theorem 2.11. Let s,t > and < 9 < /i. Then the limit in (|2.19p exists for ^ > and 

is given by 

(2.20) 

^ max {sMeiO " tM^^e{-0, tM^-eiO -sMg{-0}, < ^ < ^ A (/i - 0) 
oo, £.>0Ain-0)- 



A9,(.,t)(e) 



Remark 2.12. Let the parameters < 9 < fj, he given. The characteristic direction is the 
choice 

(2.21) (s,t) = c(^i(^ - 6'),^'i(6')) for a constant c> 0. 

With this choice the variance of logZ^^^^j ^^^^ is of order n^/^, while in other directions 

the fluctuations of log Z^^^ ^^^^ have order of magnitude n^/^ and they are asymptotically 
Gaussian |31] . By this token, we would expect the large deviations in the characteristic 
situation to be unusual, while in the off-characteristic directions we would expect the more 
typical large deviations of order e"" in both tails. In Lemma [4.2( b) we give a bound on the 
left tail that indicates superexponential decay under (I2.2ip . This also implies that if ()2.2ip 
holds, then formula (I2.20|) can be complemented with the case Ae,(s,t)iO = P^^\s,t)S^ for 
^ < 0. Presently we do not have further information about these large deviations. 

Remark 2.13. If the two sums in (|2.16p were independent we would have ^e,(s,t)iO — 
sMg{^) + tMf^-g{^). Obviously ()2.20p reflects the strong negative correlation of these sums, 
but currently we do not have a good explanation (besides the proof!) for the formula that 
arises. 

The maximum in (j2.20p comes from the choice of the first step of the path: either 
horizontal or vertical. Corresponding to this choice, define partition functions 

[nsj . k X 
k=l ^i=l ^ 
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(2-23) <irM=Efn^o.>° 



lns\,lnt\ - M i "O.J '^il/),{lns\,lnt\)^ 

£=1 ^j=l ^ 



together with l.m.g.f.'s 



7(9), hor 



(2.24) A^Zt){0= hm n-MogEe^^°*^^L-J.L"'J 
and 

(2.25) mAO= lim n-ilogEe^^°^^ 



{e),vcr 
[nsj , [ntj 



Then ^[^J.LntJ = ^KntJ + ^KntJ ^""^"^^ *° 

(2-26) A,,(,,,)(0 = A^°(^,,)(C)VA™^,)(0 

which is the starting point for the proof of ()2.20p . 

The horizontal and vertical partition functions are in some sense between the stationary 
one and the one from ()2.3p with i.i.d. weights. It turns out that these intermediate partition 
functions behave either like the stationary one or like the i.i.d. one, with a sharp transition 
in between, and this holds both at the level of the limiting free energy density and the 
l.m.g.f. Let us focus on the horizontal case, the vertical case being the same after the swap 
sot and 6* o — 6*. 

Qualitatively, with t fixed, when s is large ■^[^^jj'™ behaves like Z^^^ ^^^^ , and when s is 
small zf^^'l^™^, behaves like Z|„„i i„-n from (12.31). Here are the conditions for the transitions: 



lnsJ,LniJ oenavesuKe^Ln^J.LntJ 

(2.27) s^i(0) > t^i(^-0) 
and 

(2.28) s{^o{e) - ^o{e - 0) > i(^o(/i -0 + 0- ^o(m - ^))- 

By the concavity of ^'o and the fact that ^-i = ^'q, (fT??!) implies ([2:28]) for all ^ > 0. 
Assuming the limit exists for the moment, define 

(2.29) p(^)^^-(.,t)=Jimn-MogZfX,j. 

In this next theorem the functions p{s,t) and As^tiC) are the ones defined by ()2.5p and 

Theorem 2.14. Let s,t > 0, < 9 < n and < ^ < 0. 
(a) The limit in ()2.29p exists and is given by 

(2 SO) «W't^-rs t) = /^^'^(''*)' ^ ^"^'^^ 

^ ■ ^ ^ I ' ^ \p^s,t), z/ (EST]) 



(b) The limit in ()2.24p exists and is given by 

'sMg{0-tM^^ 

As AO, if fails. 



(2,31) ^hor ^^(^) _ J sMeiO - tM^^e{-i), if holds 
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Remark 2.15. We saw in ()2.5p that the hmitmg free energy p{s,t) of the i.i.d. model is the 
minimal free energy of the stationary models with the same bulk parameter /i. This link 
does not extend to the l.m.g.f.'s: for < ^ < As^tiO < ^9,(s,t)(C) for all 9 G (0,/x). We 
observe this at the end of the proof of Theorem 12.111 in Section \5\ 

3. The right tail rate function in the general case 

The proofs of the results for the log-gamma model utilize regularity properties of the 
rate function J of (jl.4p . These properties can be proved in some degree of generality, and 
we do so in this section. So now we consider 

(3.1) ^u= X] e^'='"^^^'^ 

as defined in the Introduction, with u E Z'^, general d > 2, and general i.i.d. weights 
{^(u)}. 

We assume 

(3.2) 3^ > such that E(e^l'^(")') < oo. 

This guarantees the existence of a Cramer large deviation rate function defined by 

n 

(3.3) I(r) = — lim lim log P-^ n^^ \^ <^(uj) € (r — e,r + e)>. 

i=l 

(Above {uj} are any distinct lattice points.) We state first the existence theorem for the 
limiting point-to-point free energy density. We omit the proof because similar superadditive 
and approximation arguments appear elsewhere in our paper, and refer to [15] . Let us also 
point out that assumption (j3.2p is unnecessarily strong for this existence result but our 
objective here is not to optimize on this point. 

Theorem 3.1. Assume ()3.2p . There exists an event fio C o/ full F-probability on which 
the convergence 

(3.4) p(y) = lim n"^ log Z^^yj 

happens simultaneously for all y S M^J.. Limit (13. 4p holds also in L^{¥). As a function of 
y, p is concave and continuous on R^j.. 

Next the right-tail LDP. To avoid issues of vanishing probabilities and infinite values of 
the rate, we make this further assumption: 

(3.5) Vr < oo : P{w(0) > r} > 0. 

Theorem 3.2. Assume ([32]) and ([33]). Then for u € M^|. \ {0} and r € M this R+-valued 
limit exists: 

(3.6) >^u('') = — lim 'T'"^ logP{log Z|„u| > nr}. 

As a function of (u, r), J is convex and continuous on (M'|_ \ {0}) x M. Ju{f) = iff 
r < p(u). 
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Let US also remark that the weight u){0) at the origin is immaterial: the limit is the same 
for so for u G \ {0} and r € M, 

(3.7) JJr) = - lim n'MogPllog ZR„ , > nr}. 

We observe this at the end of the proof of Theorem 13.21 

With a further assumption on the Cramer rate function of the weight distribution defined 
in (j3.3p we can extend the continuity of Ju to u = 0: 



(3.8) Ooo = lim X ^I{x) < oo. 



p.Sp is equivalent to requiring that I{x) < oo for all large enough x, so of course (j3.8p 
requires (j3.5p . The constant Ooo is the limiting slope of / at oo which exists by convexity. 
When assumption ()3.8p is in force we define 




(3.9) Jo(r) = 

Theorem 3.3. Under assumptions p.2p and ()3.8p . and with Jq defined by (j3.9p . Ju(?") is 
finite and continuous on x M. 

Remark 3.4. Assumption (j3.8p is in particular valid for the log-gamma model. For ~ 
Gamma(/i) the Cramer rate function for uj = logY is 

(3.10) I^{r) = -r%\-r)-logT{%\-r))+iir + logT{n), r G M. 

The limiting slope on the right is Ooo = while the limiting slope on the left would be 
limr_>._oo I'i''') = —CO. In this case Joir) is also the "rate function" for the single weight at 
the origin: 

(3.11) Jo(r) = - lim n^^logPilogF > nr}. 

ra— )-oo 



The remainder of this section proves Theorems l3.2l and l3.3| and then we prove two further 
lemmas for later use. 

Proof of Theorem \3.2l For m,n £ M+ let x^,™ € {0,1}'^ so that [(m + n)uj = [muj + 
[nuj + Xm,n ■ By superadditivity, independence and shift invariance 

IP{log^[(m+n)uJ >im + n)r} 
^ ' ^ > F{log ^L^uJ > mr}P{log > nr }P{log > 0}. 

By assumption ()3.5p there is a uniform lower bound P{log .^xm,„ > 0} > /> > 0. Thus 
t(n) = log P{log .Z^^^uj > nr} is superadditive with a small uniformly bounded correction. 
Assumption (|3.5p implies that t{n) > —oo for all n > no. Consequently by superadditivity 
the rate function 

(3.13) Ju(r) = - lim n"MogP{logZ|„u| > nr} 

exists for u = (ui, . . . ,Ud) G and r G M. The limit in (I3.13P holds also as n — )• oo 
through real values, not just integers. 



12 



N. GEORGIOU AND T. SEPPALAINEN 



Similarly we get convexity of J in (u, r). Let A € (0, 1) and assume (u, r) = A(ui,ri) + 
(1 - A)(u2,r2). Then 

n-ilogP{logZL„uJ > nr} > A(An)-i logP{log Z^^au^J > nXn} 

+ (1 - A)((l - A)n)-i logP{log %(i_A)u2j > n{l - X)r2} + o(l) 
and letting n — )■ oo gives 

(3.14) Ju{r) < XJuAn) + (1 - A)Ju2(r-2). 

Finiteness of J follows from ()3.5p . so now we know J to be a finite, convex function on 
(M^|_ \ {0}) X M. This implies that J is continuous in the interior of (M^|_ \ {0}) x M and 
upper semicontinuous on the whole set (M!J[. \ {0}) x R [28l Thm. 10.1 and 10.2]. 

The law of large numbers for the free energy implies Ju('') = for r < p{u) and then by 
continuity for r < p{u). With a minor adaptation of [9*, Prop. 3.1(b)] we get a concentration 
inequality: given u, for e > there exists a constant c > such that 

(3.15) P{|logZL„uJ -ElogZ^^ujI > ne} < 2exp{-ce'^n) for ah n e N. 

Since n^^Elog Z|^„uj — >■ p{u), this implies that Juir) > for r > p{u). 

We do a coupling proof for lower semicontinuity. Let (u, r) — )• (v, s) in (M^J. \ {0}) x M. 
If each coordinate Vi > then we have continuity J-air) — )■ Jv{s) because convexity already 
gives continuity in the interior. Thus we may assume that some coordinates of v are 
zero. Since coordinates can be permuted without changing J, let us assume that v = 
{vi,V2, ■ ■ ■ ,Vk,0, ■ ■ ■ ,0) for a fixed 1 < k < d where vi, . . . ,Vk > 0. If eventually u is also of 
the form u = (ui,U2, . . . , Ufc, 0, . . . , 0) for the same k then we are done by convexity-implied 
continuity again, this time in the interior of (M^ \ {0}) x M. 

The remaining case is the one where ui, . . . ,Uk > and (ufc+i, . . . , Ud) 0. We develop 
a family of couplings that eliminates these d — k last coordinates one by one, starting 
with Ud, and puts us back in the interior case with continuity. Denote a lower-dimensional 
projection by ui^^ = (ni,n2, . . . ,'Ufc). 

The set of paths '^\nvi\ is decomposed according to the locations of the [nti^J unit jumps 
in the e^-direction. The projections of these locations form a vector vr from the set 

ALnuJ = {vr = {x*}£o'^+' € (4-i)L-.J+2 : o = < < • • • < xL^'^J+i = Lnui,,_ij}. 

The partition function then decomposes according to these jump locations: 

(3-16) Z\nu\ = X] ^(0,0),{xi,0) n -^(x',i),(x«+i,i) = X] 

where the last equality defines the d — 1-dimensional partition functions Z^^. 
For a fixed vr, define a new environment uj indexed by with this recipe: 

(i) For < i < [nn^J: for y e Z'fr^ such that x* < y < x*+i but y 7^ x\ set 
'^(y) = w(y,i). 

(ii) a;(0) =a;(0,0) and for 1 < i < \nud\, uj{\nu.i^d-i\ + i^d-i) =a;(x*,i). 

(iii) Pick all other tj(y) independently of everything else. 



(3.17) 
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Now, keeping tt fixed, we project the paths down to and create a partition function 
(marked by a tilde) in the new environment u). 

[nuci\ 

logZ^ = logZ(o,o),(xi,o) + Yl log^(x%i),(x«+\i) 

i=l 

= Y logZ(x»^i)_(x»+i,i) + ^ w(x^^i) 

i=0 i=l 

= Y logZ^^^x'+i + ^ a;([nui^rf_ij + ied-i) 

< log^[„ui,d_iJ + Lnnded_iJ- 

Introduce the continuous functions [1 <i < d) 

i-l 

(3.18) Fj(u) = ^ [{uj + Uj) log(nj + Ui) - uj logUj - Ui logUj). 

i=i 

Counting the number of ways to decompose the length from to \nui\ into [nn^J + 1 
segments and Stirling's formula give 

mo = \\nu^= n f^T^^tT!"'^) =exp{nF,(u) + o(n)} 

(3.19) i<fiiV + 1 y 

< exp{nFrf(u) + nr]} 

where the last inequality is valid for large n and we introduced a small r/ > that we can 
send to zero after limits in n have been taken. By a union bound and the coupling (j3.17p 
separately for each tt G A|^„uJ , 



-J\iif) < lim n ^ log P| log > nr — log mo} 



TreA 



[nuj 



-ni^(^^T^ + '' ^logp|logZL„ui,,_iJ + Ln«,e,_J > nr - nFd(u) - n7?|^ 
= Fd(u) - Ju^^_j+„^ed_i(r-Frf(u) -r/). 

In the last step above a little correction as in (|3.12|) replaces [nui + \nU(iGd-i\ with 
[nui^d-i + nuded-i\. 

Let Mi^d = u and for 1 < i < d, 

d 

j=i+i 

Proceeding inductively, we get the lower bound 

(3.20) Ju(r)> Js,,,(r- {F^i^) - v)) - ^ ^^(")- 

k+l<i<d k+l<i<d 
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On the right-hand side we have a rate function Jq^ ^ with ui_fc — > vi^/j in the interior of 
M'^. Thus we have continuity. We can first let rj\0. Then let (u, r) — )• (v,s). Note that 
tij — 7- implies Fi{u) — )• 0. Together all this gives the lower semicontinuity: 

hm Juir) > Jv^,^:is) = Jv{s). 

(u,r)— >(v,s) 

Now we know J is continuous on all of (M^j. \ {0}) x M. 

Let us observe limit (j3.7p . From one side we have 

P{logZp„„j > nr} > P{logZL„uJ > nr}P{w(0) > 0}. 

From the other, pick a coordinate Ui > 0, and for each n an integer n < nin < n + o{n) 
such that 2ej + [nuj < [m^uj. For each n fix a directed path {x^} from 2ei + [nuj to 
[m„uj. Inequality 

i 

gives 

P{logZ° ^j > nr}F[u{e,) + J]c.(x,") > o} < P{log > nr} 

j 

Assumption (j3.5p and the continuity of J give the conclusion. □ 

Proof of Theorem VJ.'Ji It remains to prove continuity at (0,s). Let (u, r) — )• (0,s). Define 
the right-tail Cramer rate function for a > 0, x € M: 

[naj 

Ka(x) = — lim n^^ logWl n^^ wfxj) > nx> 

i=l 

Check that as {a,x) — )• (0, s), Ka{x) Jo(s) defined by <\3.9^ . 
For upper semicontinuity, bound below by a single path: 

ILnuJIi 

Ju{f) < — lim n^"*^ logP< n^^ > <^(xi) > nr > = Ki^\Ar). 

i=l 

For lower semicontinuity, permute the coordinates so that ui > as u ^ 0. Apply ()3.20p 
after r] has been taken to zero: 

Ju(r)> J„,ei(r- J2 ^*("))- J2 ^*(")- 

2<i<d 2<i<d 

Since Ju^ei = we get the lower semicontinuity. □ 

Finally two lemmas for later use. The next one allows more general lattice sequences for 
the right-tail LDP. 

Lemma 3.5. Let y € (0,00)'^ and u„ € he a sequence such that n^'^Un y- Then for 
r G M, 

(3.21) lim n-MogP{logZu„ > nr] = -Jy{r). 



al{x/a), X > aE[w(0)] 
0, x< aE[w(0)]. 
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Proof. Let us use assumption (I3.5p again. Since the coordinates of u„ and [nyj are increas- 
ing to CO, for each n we can find in and such that [^nyj < u„ < [m„yj and in such a 
way that n — £n, n — nin are eventuahy o(n). For each n fix directed paths {xn,i}Q<i<Kn 
from [lny\ to u„ and {x'nj}o<j<K'^ from u„ to [m„yj. Then 

where 

logW„= ^ uj{xi) and logl<= ^ a;(2;-). 

l<i<Kn ^<i<K'^ 

Assumption n^^u„ ^ y imphes that Kn and K'^^ are also o(n). 
The estimates we need fohow. For example, 

P{logZL„„yj > nr} > P{logT< > 0}P{logZu„ > nr} 

and then by assumption ()3.5p and the continuity of the rate function, 

lim logP{log > nr} < lim logPilog Zi^^ yi > nr} = —Jy(r). 

Similarly for the complementary lower bound on lim. □ 

Lemma 3.6. Suppose that for each n, Ln and Zn are independent random variables. Assume 
that the limits 



(3.22) A(s) = - lim n"^ logP{L„ > ns}, 

(3.23) cl){s) = - lim n"^ logP{Z„ > ns} 

exist and are finite for all s G M. Assume that X{ax) = (j){a(f,) = for some ax,a(j) € M. 
Assume also that X is continuous. Then for r G M 

(- inf {(/)(r - s) + A(s)}, r > as + a\ 

(3.24) lim n~^ logP{L„ + Z„ > nr} = { «A<^'<»--a^ 

(O, r<a^ + ax. 

Proof. The lower bound > follows from 

P{L„ + Z„ > nr} > P{L„ > ns}P{Z„ > n{r - s)}. 

Since an upper bound is obvious, it remains to show the upper bound for the case 
r > + ax. Take a finite partition ax = qo < • • • < Qm = r — a^. Then use a union bound 
and independence: 

P{i„ + Zn> nr} 

< F{Ln + Zn> nr, Ln < nqo} 

m— 1 

+ ^ P{L„ + Zn> nr, nqt < Ln < nqi+i} + P{L„ > nqm} 

1=0 

m—l 

< F{Zn > n{r - go)} + Yl ^ "(^ " %+i)}IP{^n > nqi} + P{L„ > nq^}. 



i=0 
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From this 

lim n^^ logP{L„ + Z„ > nr} 

n— >oo 

< -m.\n\<f>{r - qo), min [(f){r - qi+i) + X{qi)], X{qm)\ ■ 

K 0<i<m— 1 J 

Note that \{qo) = (j){r — q^n) = 0, refine the partition and use the continuity of A. □ 
4. Proofs for the i.i.d. log-gamma model 

In this section we prove the results of Section [2. II Throughout this section the dimension 
d = 2 and the weights satisfy Y-~^ ~ Gamma(/i) as in (|'2.1|) . As before, for {s,t) G 
\ {(0,0)} define the function Jg^t by the limit 

(4.1) Js,t{r) = - lim n"! logP{log ZL„,j,LniJ > nr}, r G R. 
At the origin set 

(4.2) A.M = ("' 

I fir, r > 0. 

Then, as observed in Remark l3.4l the function Js,tif) is finite and continuous at all (s, t, r) S 
X M. 

We begin with a lemma that proves Theorem 12.11 
Lemma 4.1. For {s,t) G the limiting free energy of (12. Sp satisfies 

(4.3) p{s, t) = inf {-s^oiO) - t^o{^^ - 9)}. 

o<e<fi 

The infimum is achieved at some 9 because ^'o(O-l-) = — oo. 

Proof. The proof anticipates some themes of the later LDP proof but in a simpler context. 
We already recorded the law of large numbers (I2.18P . The decomposition (see Figure [2]) 

\ns\ . k V \nt\ . t . 

(4-4) ^M.LniJ = 5Z ( n^^.O ( n^O'i )^{°Q,(M,L"^J)■ 

from (|2.14p gives asymptotically 



n 

{fc,i),{L"^J,L"tJ) 



^lim log = „1™o{i<^|l-J |J ^ ^' 

V l<T|Ln*J ' ?5 ^ " ' ^Pl/),(LnsJ,LntJ)) } • 

This can be coarse-grained with readily controllable errors of sums of independent variables. 
We omit the details since similar arguments appear elsewhere in the paper. The conclusion 
is the alternative formula 

Jim^ n log ^ ^^^j 

(4-5) , , 

= sup {-a^'o(6')-Fp(s-a,t)}\/ sup {-6^'o(^ - 6^) p(s, t - 6)}. 

0<a<s 0<6<i 



LARGE DEVIATIONS FOR LOG-GAMMA POLYMER 



17 



[nt] 




1 



[nsj 



Figure 2. Graphical representation of the decomposition in equation ()4.4p . 

Take s = t, combine (|2.18p and (|4.5p . and use the symmetry p{s, t) = p(t, s) to get 
-t{^o{e) + ^o{fi-e)) = sup {-a(^'oWA^o(/i-^))+p(t-a,t)}. 

0<a<t 

Take 6 E (0,/i/2] so that ^'o(^') < 'J'o(/^ — ^) [^o is strictly increasing] and set a = t — s: 

-mo{fi-e)= sup {s^o{9)+p{s,t)}. 

0<s<t 

Turn this into a convex duality through the change of variable v = ^q{9): 
(4.6) -t^o{t^-%\v)) = sup {sv + p{s,t)}, € (-oo,^'oW2)]. 

0<s<t 

It follows from the limit definition of t) that it is concave and continuous in s € [0, t]. 
Extend /(s) = —p{s,t) to a lower semicontinuous convex function of s G M by setting 
f{s) = oo for s ^ [0,t]. Then gSD tells us that 

nv) = -t^o{^^-nHv)) for (-oo,vI/oW2)]. 

From this we can differentiate to get lim„\^_oo(/*)'(^^) = and q{h/2)) = t. These 

derivative values imply that for s € [0, t] the supremum in the double convex duality can 
be restricted as follows: 

f{s)= sup {vs-f*{v)}. 

i'G(-oo,*o(At/2)] 

Undoing the change of variables turns this equation into (j4.3p which is thereby proved. □ 

The next lemma gives left tail bounds strong enough to imply Is,t{f) = oo for r < p{s, t), 
and the same result for the stationary model. The proof is a straightforward coarse-graining 
argument. We do not expect the results to be optimal. 
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Figure 3. The [msj x [mt\ rectangles and the diagonals Aj in the proof 
of Lemma 14.21 The thickset line is a lattice path that is counted in Zi . 

Lemma 4.2. Fix < a < 1. Then there exist constants < c, C < oo that depend on the 
parameters given below, so that the following estimates hold. 

(a) For (s, t) € (0, oo)^ and r < p{s, t), 

(4.7) P{logZL„,j,LntJ < nr} < Ce"^"'^" for all n>l. 

(b) For {s,t) = Q(^i(/i — for some a > 0, parallel to the characteristic 
direction, and r < p^^\s,t), 

(4.8) P{log zf2j,Lntj < < Ce-^"'^" for all n>l. 

Proof. We give a proof of (b) with some details left sketchy. Part (a) has a similar proof. 

(6) 

We bound ^^t] below by considering a subset of lattice paths, arranged in a 

collection of i.i.d. partition functions over subsets of the rectangle. 

The choice of (s,t) implies that p^^\s,t) = p{s,t). Fix < e < {p^^\s,t) - r)/4. Fix 
m € N large enough so that m{s At) > 1 and 

(4.9) ElogZL^^j^L^tj > m(r + 2e). 

Let i?^'^ = {a, . . . ,a + k — 1} x {b, ...,£ + b — 1} denote the k x £ rectangle with lower 
left corner at {a,b). For i,i >0 define pairwise disjoint [msJ x \rnt\ rectangles 

Define a diagonal union of these rectangles by Aj = U^>o i > 0. (See Figure [3|). 

Let M = [n^\ [ms\ . This is the range of diagonals Aj we consider. Then we cut the 
diagonals off before they exit the [nsj x \nt\ rectangle. Let N = N(n) be the maximal 
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integer such that Hes in [0, \ns\] x [0, \ nt\]. Diagonal Aa/ exits the \ns\ x \ nt\ rectangle 
through the east edge, and consequently there exist positive constants Cm, Cm such that 

(4.10) [nsj -Cm< N[ms\ + [ms\ [n"J < [ns\, and [nt\ - Cmn" < N[mt\ < [nt\. 

Having defined the cutoff N, define the other diagonals by A" = IJo<£<Af ^ — 

i < M. These diagonals lie in [0, [ns\] x [0, [nt\]. Fix a path vr that proceeds horizontally 
from {N[ms\,N[mt\ + 1) to ([nsJ, A^[mtJ + 1) and then vertically up to ([nsJ, [?T.tJ). The 
number of lattice points on vr is a constant multiple of n'^. 

For < i < M let Zj denote the partition function of paths x. of the following type: 
rc. proceeds along the x-axis from the origin to (i[msj + 1,0), enters A" at (i[msj + 1, 1), 
and stays in A" until it exits from the upper right corner of B^j^ with a vertical step that 
connects it with tt. After that x, follows vr to ([rasj, [ntj). The number K of points on x, 
outside A" is independent of i and bounded by a constant multiple of n". Let 

X = min{ya; : X € vr or X G {{i, 0) : < i < M}} 

be the minimal weight outside Af encountered by any path x. of Zi, for any < i < M. 

Let be the partition function of all lattice paths in A" from the lower left corner of 
Bq to the upper right corner of i?]^- Then Zi > Zf^, and consequently 

M 

P{logZ[2j,L„ij < ^r} < P{ log^X^Zf < nr} 

M 

= pjxiogX + log^Z^ < nrj 
i=0 

M 

(4.11) <¥{K\ogX < -ne} + p| log ^ < n(r + e)|. 

i=0 

Explicit computation with the gamma distribution and K < en"" give T{K log X < —ne} < 

2 

e~" for large n. 

The {Z^} are i.i.d., and Z^ is a product of the i.i.d. partition functions Z^ of the 
individual rectangles B^ whose mean was controlled by ()4.9p . A standard large deviation 
estimate for an i.i.d. sum gives 

pjlog^Zf <n(r + e)| <P{logZ^ <n(r + e)}^'^ =p|^logZ^ <n(r + e)y 

1=0 k=0 

n/m+o{n) 

= P{ log Z^<nir + e)y < e"""^'^ < £"'=1"'''". 

k=0 

Putting these bounds back on line (j4.1ip completes the proof of (j4.8p . □ 

The main work resides in proving the following right tail result. 
Proposition 4.3. Let {s,t) £ M^. Then for all r € M, Js,tir) is given by 

(4.12) JsAr)= sup {re- inf (tMg{0 - sM^_e{-0)} ■ 



M 
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Before turning to the proof of Proposition 14.31 let us observe how Theorem 12.21 fohows. 

Proof of Theorem \2. 21 Only a few simple observations are required. Start by defining Ig^t 
as given in ()2.6p . Then formula (j2.9p that connects Is,t and Js^t is established by (j4.12p and 
by knowing that Js,tii^) = for r < p{s,t) (Theorem I3.2p . The regularity properties of Ig^t 
follow from the general properties of J in Theorems 13.21 and 13.31 

The upper large deviation bound (j2.8p is built into (j4.7p and (j4.ip . 

For the lower large deviation bound (12. 7p we consider three cases. 

(i) If pis, t) G G then P{n-Mog Z^n^j^LntJ G G} ^ 1 and ([231) holds trivially because its 
right-hand side is < 0. 

(ii) If G C (—oo, p{s,t)), (12.7j) holds trivially because its right-hand side is —oo. 

(iii) The remaining case is the one where G contains an interval (a, 6) C {p{s,t), oo). 
Since the distribution is continuous including a into G makes no difference, and so 

n~ilogP{n-MogZL„,j,LntJ G G} 

> n^^log (P{logZL„sj,LntJ > na} - ¥{log Zins},[nti > nb}) 
— -JsA"-) 

where the limit follows from (14. ip and the strict increasingness of Js^t on \p{s,t), oo) which 
implies that for large enough n 

P{logZL„,j,L„tj >n6}<e-"^-*('^)-"^ 

for some e > 0. We can take a = inf G fl t), oo) and then Jg j(a) = inf Is t(r) 

' reGn{p(s,t),oo) ' 

The remainder of the section is devoted to proving Proposition 14.31 Again we begin with 
the decomposition (I4.4p of the stationary partition function. Inside the sums on the right 
of (|4.4p we have partition functions with i.i.d. Gamma~^(/i)-weights {i^j} whose large 

deviations we wish to extract. But we do not know the large deviations of logZ^^^^j ^^^^ 
so at first the decomposition seems unhelpful. To get around the problem, use definition 
(|2T5]1 to write 

[nsJ 

log ^\nli , [nt] - ^Ant] = Yl L^tJ " 

i=i 

By Proposition 12.10] we have a sum of i.i.d. 's on the right, whose large deviations we can 
immediately write down by Cramer's theorem. To take advantage of this, divide through 
(Hai) by Z^^^l^^ = Yl\=l Voj to rewrite it as 

[nsJ [ntj , [nij 

/),(LnsJ,LntJ) 



(4.13) 



1=1 e=i ^ j=e+i ^ 

[nsJ , [ntJ \ / \ 

+ { n ^07/ ) ( n ^^'0 ) ^g,l),([nsJ,LniJ)- 
fe=l ^ 7 = 1 ^ ^ i=l ^ 
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To compactify notation we use a convention where the y-axis is labeled by negative indices 
and introduce these quantities: 



(4.14) 



for A; € Z, % = < 



llj=-k+l ^0,j ' 



^ ( ni=J v^o^ ) nil u,o, 

where an empty product equals 1 by definition, and 



k<0, 



A; > 1, 



for z € 



v(z) = < 



(4.15) 

Then (|4.13p rewrites as 
(4.16) 

from which we extract these inequalities: 



{l,l-z\), z<-l, 
(1,1), -l<z<l, 
[{[z],l), Z>1. 



[nsj [nsj 
i=l k=-lnt] 



(4.17) 



log^fc + log^v(fc),(LnsJ,LntJ) 



[ns] 

< ^logC/j^LniJ 



i=l 



These inequalities will be the basis for proving Proposition 14.31 

We record the right tail rate functions for the random variables in ()4.17p . 

For the i.i.d. weights {C/j,[ntj} we have the right branch of the Cramer rate function 



(4.18) Rs{r) = - Jim^n-ilogpj |^log[/,,LntJ > nrj = 



r < -s^'o(6'). 



The rate function Ig defined by (j4.18p is given by 

(4.19) /e(r) = -r^'o^(-r) -logr(^'o ^(-r)) + 6lr + log r(6'), r E M. 

The convex dual of Rs is given by 



(4.20) 



siogr(0-O -siogr(6'), o<^<9 

oo, ^ < or ^ > 



and we emphasize that it can be finite only when 9 > > 0. 
For real a G [—t, s] 



(4.21) 



Knir 



lim n ^ logP{log r/|„„| > nr} 
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exists and is finite, convex and continuous in r. (For a < it is simply a Cramer rate 
function for an i.i.d. sum, and for a > we can use Lemma 13.61 ) The convex dual is 

<(6 = sup{^r - Ka{r)} 

'(t + a)(logr(^-0 + e) -logr(/u-0)), -t<a<0, ^ > 0, 

t ( log T{fi-e + o- log r(^ - e)) + a( log T{e - - log r{e)) 

0<a<s, 0<^<0, 

^oo, otherwise. 

The derivation of (I4.22p is similar to that of (I4.20p from (I4.18p . Note that there is a 
discontinuity in and k* as a passes through 0. The rightmost zero m^^a of Ka is the law 
of large numbers limit: 

logr?[naJ _ /(t + a)^'o(/U-6l), -t<a<0 



(4.22) 



(4.23) m«;,a = lim 



n I t^o(/i - 6*) - a^o(6'), < a < s. 



In contrast to the functions Ka and k*, m^^a is continuous at a = 0. Introduce the "macro- 
scopic" version of (|4.15p : for real a. 



(4.24) n V(na) v(a) 



(0, -o), -t<a<0, 
(a,0), 0<a<s. 



With this notation we have, again for real a G [— t, s], for the partition functions that appear 
in (j4.16p these large deviations: 

(4.25) J(,,i)_v{a)(?-) = -^lim^n-^logP{logZ°(„„)^(L„,j_L„^j) > nr}. 

We used Lemma [331 to take care of the small discrepancy between ([nsj, \ nt\) —w{na) and 
[n((s,t) — v(a))J, unless a = —t ov a = s when this is a case of i.i.d. large deviations, and 
therefore simpler. 

Let mK,a and mj^b be the rightmost zeroes of Ka and J(s,t)-v(b) respectively. For (a, b) e 
[-t,s]2 let 

Hs^ir) = lim n-^logP{log??LnaJ + log ^v(n6),{LnsJ,LntJ) > ^r} 



(4.26) 



r < m^^a + mj^b 

inf {Ka{x) + J(s,t)_v(b)('' - x)}, r > m^^a + mj^b- 



The existence of H'^'^{r) and the second equality follow from Lemma 13.61 We need some 
regularity: 

Lemma 4.4. Fix < s,t < oo and a compact set K C Then Hg'^{r) is uniformly 
continuous as a function of (6, r) G [—t, s] x K , uniformly in a & [—t, s]. That is, 

(4.27) hm sup \H:f{r)-H:f{x)\ = 0. 

a,b,b'e[-t,s],r,x&K: 
\b-b'\<S ,\r~x\<S 
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Proof. This follows from the explicit formula in (14.260 . First, we have the joint continuity 
(6, r) I—)- ()_v(6)(?^) from Theorem 13. 31 Second, we argue that x in the infimum can be 
restricted to a single compact set simultaneously for {a,b,r) £ [—t,s]'^ x K. That m^.a is 
bounded is evident from (|4.23p . To show that the upper bound r — mjf, of x is bounded 
above, we need to show a lower bound on mjb = p{(s,t) — v(5)). A lower bound on the 
free energy is easy: by discarding all but a single path, 

p((s,t)-v(6)) =^lim n-ilogZ°((^_,)_-(,))j > -{s + t - \b\)^o{f^). □ 
We abbreviate H^tir) = H^'^tir). 

The unknown rate functions Js^t are now inside (j4.26p . while the other rates Rg and Ka 
we know explicitly. The next lemma is the counterpart of ()4.17p in terms of rate functions. 

Lemma 4.5. Let s,t>0 and r G M. Then 

(4.28) Rsir) = inf Hl.ir). 



Proof. For any a € [—t, s], by the first inequality of (j4.17p . 

-i?s(r) = lim n~MogP<^ log f/j i ntj >nr\ 

i=l ■' 
- i™o""^^°S^{ log^L^aJ +log^v(na),(LnsJ,LntJ) > 

(4.29) > -Hl^ir). 

Supremum over a € [—t, s] on the right gives < in (I4.28|) . 

To get > in (I4.28P we use the second inequality of (I4.17P together with a partitioning 
argument. Let e > 0. Note this technical point about handling the errors of the partitioning. 
With B,5 > 0, Chebyshev's inequality and the l.m.g.f. of (j2.2p give the bound 

(4.30) p{X^logy,i < -ne\ < e""^!^"^"^^!"^ ^) < g-^^/^ 

^ i=i ^ 

where the second inequality comes from choosing 5 = 6{e, B) small enough. The right tail 
for logy does not give such a bound with an arbitrarily large B. Consequently we arrange 
the errors so that they can be bounded as above. 

Given i? > 0, fix a small enough 5 > and let — t = ao < ai < • • • < = < • • • < 
a^a = s be a partition of the interval — [t,s] so that |aj+i — aj| < 5. We illustrate how a 
term with index k from the right-hand side of (|4.17p is reduced to a term involving only 
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partition points. Consider the case > and let [najj < k < [na^+ij. 

P{ logr/fc + log ^v(fc),(LnsJ,LntJ) > 



, lnai+i\ k-1 

< P| logr/L„a,+iJ + log^°(„„^),(L„,j,LntJ) - - ^ logl"j,i > nr [ 

< P{ logr/[na,+ iJ + log^v(na,),{Ln^J,M) ^ " ^)} 

+ Pi - ^ logC/j,o- logyi,i>nel 
^ 1, 1 1 I — I ' 



j=k+l j=L"aiJ 
(4.31) < F{ log7?Lna,+,J + log^°(„,,),(L„.J,LntJ) > ^(^ " ^)} + « 



-Ben/2 



On the other hand, if Oj < and [— noj+ij < —k < [— nojj, then we would develop as 
follows: 

log^fc+log^v(fc),{Ln.J,LntJ) 

-Inai] -k-1 

<logr?Lnad- logFoj +log^v(na,+ i),{Ln.J,LntJ) - Yl ^Og^l-i 

j=-fc+l j=[-na,+iJVl 

and get the same bound as on line (I4.3ip but with Oj and Oj+i switched around. 

Now for > in (j4.28p . Assume n is large enough so that ne > log(ns + nt). Starting from 

(BIZl) 

lns\ 

n ^iogr<; "S^ loe:Uii„.t\ >nr 



-ilogpj J^logC/i,L„tJ > 
1=1 



- , .^i^^i ^logP{log^fc + log^v(fc),(LnsJ,LntJ) >^('^-^)}+'^ ^ \og{ns + ut) 

— \nt\<k<\ns\ \ jiL J/ 

- 0^"'^ (^^ ^"^''L-d + log^v(na,+o,(Ln«J,Ln*J) ^ ^^^^ " 2^)) + e"^^"/') 

Take n — >■ oo above to obtain 

-Rs{r) < { ^ nmx_^ ( - F^f (r - 2e)) V ( - Be/2) } 



(jr<j<r7i— 1 

< sup ( - H'^fir - 2e)) V ( - Be/2) + e. 

a,be[-t,s]: \a-b\<S 

We first let J \, 0, and by Lemma 14.41 the bound above becomes 

-Rsir) < sup ( - H'^'tir - 2e)) V ( - Be/2) + e. 

a&[—t,s] 
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Next we take B oo, and finally e \ with another application of Lemma 14.41 This 
estabhshes > in (i428]l . □ 

A key analytic trick will be to look at the dual <^(* f)_v(a)(^) right tail rate as a 

function of a. This lemma will be helpful. 



Lemma 4.6. For a fixed ^ € [0, fi) the function 

f-'^(W^(a)(0, ae[o,t] 
\ oo, a < or a > t 



(4.32) G5(a) 



is continuous on [0, t], and convex and lower semi- continuous on M. In particular, G|*(a) = 
G^{a) forae M. 

Proof. To show convexity on [0, t] let A € (0, 1) and a = Xai + (1 — A)a2: 

-•^(t,t)-v(a)(0 = -sup{er - J(t,t)-v{a){r)} 

< inf inf {X{Jt-a,An) - en) + (1 - X){Jt-a,,t{r2) - er-2)} 

reM (ri,r2): 

Ari+(1— A)r2=r 

= inf {A(Ji_,,,t(ri)-eri) + (l-A)(Jt_,„j(r2)-er2)} 
= A inf {Jt_ai,t(n) - Cn} + (1 - A) inf {Jt-a2,t(r2) - Cr2} 

(4.33) = -xj:^a,AO - (1 - A)j;_,,,t(0- 

The inequality comes from the convexity of J in the variable (t — a, t, r). 

For finiteness on [0,t] it is now enough to show that G^{a) is finite at the endpoints. 
Continuity then follows in the interior (0,t). First take a = t. Then Jq^ is the dual of a 
Cramer rate function, and for ^ > 

(4.34) G^{t) = -Jo*i(e) = -tlogEe«i°*^^^.° 

which is finite for ^ < /i. 

Convexity of Js,t{r) and symmetry Js,t{f) = Jt,s{f) imply Jt^tif) < JG,2t{f)- From this 

G'5(0) = -JM(0 = mf{JMW-M 

(4.35) < inf { Jo,2t(r) - ir} = -J^^^tiO < oo. 

rglR 

Continuity at a = 0. To show that is also continuous at the endpoints, we first obtain 
a lower bound. For any r G M, 

J:-aAO>rC-Jt-a,t{r) 
hence, by continuity of Js^t in the (s, t) argument, 

(4.36) ]mJ:-a,ti^)>rC-JtAr)- 

Supremum over r gives lim„ ^n Jt-a,tiC) ^ "^ttiO- 
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For the upper bound, let < a < t. Varadhan's theorem (Thm. 4.3.1 in [12j) apphes 
in the present setting. This is justified in the proof of Cor. 12.51 below and another similar 
justification is given for (|5.4p below. Consequently 

JttiO= lim n-MogEe«'°s^L"tJ,L"tJ 

' n— ^oo 

> lim n"MogEe^'°s^L"(t-<>)J,LntJ + Mm n^Mog Ee^^^=L"(t— )J+i '"^^-^'LntJ 

n—^oo n—^OD 

(4.37) = j;_,,,(e) + alogEy«. 

Taking a \ yields continuity at o = 0. 

Continuity at a = t. The lower bound follows as in the previous case. For the upper 
bound we use a path counting argument. Let e""^^*'*^ be an upper bound on the number of 
paths in n^„gj^^„jj such that F{0+,t) = 0. Consider first the case where < ^ < 1. Then 

lnt\ + ln{t~a)] 

iO = lim log E ( Yl n 

(4.38) [nt\ + ln{t-a)] 

< lim n-^log y TT E(y)^ 

= F(t-a,t) + (2-a/t)Jo*i(e) 
For 1 < ^ < ;U Jensen's inequality yields 

(4.39) J:_,/0 < ^F{t - a,t) + (2 - a/t).!^,^^. 

Let a t to get the continuity. 

= is a consequence of convexity and lower semicontinuity, by [28| Thm. 12.2] □ 



Proof of Proposition \4-3\ The remainder of the proof is convex analysis. The goal is to 
derive this formula for the right tail rate function Jg^t- 

(4.40) JsAr)= sup {rC- inf {tMe{0 - sM^_e{-0)} ■ 

We begin by expressing the explicitly known dual Rs{£,) from (I4.20p in terms of the 
unknown function J(s,t)-v{a)- Equation (I4.26P says that H^^ is the infimal convolution of 
Ka and J(s^t)-v{a)^ in symbols Hg^ = KaOJ(^s,t)-v{a)- By Theorem 16.4 in [28] addition is 
dual to infimal convolution. Starting with (|4.28p we have 

KiO = sup sup {r^ - (ft;anj(s,t)_v(a) )(?■)} 

-t<a<s reR 

(-4.41) = sup (KanJ(s,t)-v{a))*(0 

^ ' —t<a<s 

= sup {<{0 + Jis,t)-^ia)iO}- 

—t<a<s 

Combining this with (j4.20p gives, for < ^ < 0, 

(4.42) slogT{e-0-slogr{e)= sup {<(e) + J(Vv(a)(0}- 

—t<a<s 
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Now regard ^ G [0, /i) fixed, and let 9 E {S,, fi) vary. Introduce temporary definitions 
-hi:{e) = M^^e{-0 = log T{^l-e + i)- log T{^i-e), -t < a < o 
di:{e) = MeiO = logr(0 - - logr(0), O < a < s. 

Substitute ([02]) and (|i:i3]) into equation (^^1^ to get 

r(^) ^ r(/i-e) _t<,"<. 



(4.43) uaie) 



(4.44) .log _ t log y^/;^^^ = sup {an.(0) + J(;,)_^(,)(0} 



t(d^(0) + /ig(e)) 



The right-hand side begins to resemble a convex dual, and will allow us to solve for Jg^t- We 
can specialize to the case s = t because (t, t) — v(a) gives all the pairs (s, t) with < s < t. 
When s = t, the Js^t = -'t.s symmetry allows us to write (|4.44p as 

tid^ie) + h^ie)) = sup {a{h^i9) V d^(0)) + j;_,,,(0} 

0<a<t 

and it splits into cases as follows: 

rsupo<,<i{a/i5(0) + Jt„,t(0}, [(/i + 0/2,^) 

lsupo<a<t {a(ie(^) + JLaAO}, G (^, (^ + i)/2]. 

We can discard one of the branches above. For ii 9' = fi + — 9 then d^{9') = h^{9) and we 
see that the two equations given by the two branches are in fact equivalent. So we restrict 
to the case 9 £ [{fi + 0/2) ^) and continue with 

(4.45) tid^i9) + h^{9)) = sup {ah^{9) + J^^^AO} ■ 

0<a<t 

The function is strictly increasing so we can change variables via v = h^{9) between the 
intervals 9 € + S,)/2, n) and v € [h^{{fi + S,)/ 2), oo). Recah also G^{a) = — Ji*_at(0 from 
Lemma 14.61 This turns (j4.45p into 

(4.46) t{{d^ o h~^){v)+v)= sup {av-G^{a)} = GUv), h^{li±^) < v < oo. 

0<a<t 

Utilizing = we get an expression for the rate function J: 

(4.47) Jt^aAr)= sup K- j;_,.i(0}= sup K + G5(a)} 

= sup {r^ + sup[af — G|(w)] } 
5e[o,M) fGK 

= sup {r^ + sup [aw — G|(w)]} 

€G[0,m) i,G[/i5((m+0/2),oo) 

(4.48) = sup {r^+ sup [{a - t)v - td^{h7^ (v))]} . 

€e[0,M) i>G[he((M+0/2),oo) 

In the next to last equality above we restricted the supremum over v to the interval v G 
[h^{{lJ- + 0/2), oo). This is justified because G| is convex, a > 0, and from ()4.46p we can 

compute the right derivative {G'^)' {h^{^^^) + ) = 0. The restriction of the supremum then 
allows us to replace G|(w) with (|4.46p . 
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The proof is complete. In the case < s < t take a = t — s on line (I4.47p . Lme ()4.48p 
is the desired representation for Jg^t- It turns into (|4.40p by the v to 9 change of variable. 
The case s > t follows from the symmetry Js^tif) = Jt,si'>')- D 

The next lemma makes explicit the formula(s) for J* ^ that were implicit in the proof of 
Proposition 14.31 

Lemma 4.7. Let s,t>0 and ^ G [0,n)- Then 

(4.49) J:/0= inf {tMp(e)-sAVp(-0} 

(4.50) = ini {sMeiO-tM^-ei-O}- 



Proof. (|4.50p comes from (j4.49p by the change of variable p = p + ^ — Comparison of the 
two shows that we can assume s <t. To prove ()4.49p for s <t, start from Lemma 14.61 

J:A0 = -G^t -s) = -Gl*it -s) = - sup{(t - s)v - Gliv)}. 

Restrict the supremum as in (|4.47p - (|4.48p . substitute in (j4.46p and change variables from 
vto6 = h~'^{v). □ 

Proof of Corollary \2. 51 If ^ > ^, 

eiogz 

j 

for any particlar path x. G '^[ns\,[nt]j ^^'^ then Ag^tiO = oo comes from M^(^) = oo from 

Let ^ < /U. Pick 7 > 1 such that 7^ < p. Then the bound 

supn"MogEe'^^^°s2L'"j,Lntj < 

n 

follows from pathcounting, as in (j4.38p - ()4.39p . This bound is sufficient for Varadhan's 
theorem (Thm. 4.3.1 in [12|) which gives 

lim AgtiO = n"^logEe^'°s^L-J.LntJ = 7,%(0 = sup{r^ - t(r)} 
= sup {r^ - Ig^t{r)} = sup {r^ - Js,t{r)}. 

r>p(s,t) 

We discarded {Is,t = 00} = {r < p{s,t)} from the supremum. Since Ig^t increases for 
r > p{s,t), the case ^ < of (|2.1ip follows. For ^ > the values Js,t{f) = for r < p{s,t) 
can be put back in because they do not alter the supremum. Consequently As^tiO = Js,tiO 
for ^ > 0. Lemma 14.71 completes the proof of this corollary. □ 

There is nothing new in the proof of Corollary 12.81 so we omit it. 
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5. Proofs for the stationary log-gamma model 
In this section we prove the results of Section 12. 2i 
Proof of Theorem \2.14\ Coarse-graining and simple error bounds readily give this limit: 

= sup {-a'^oiO) +p{s - a,t)} 

0<a<s 

= sup inf {-a^o{0) + {a-s)^o{p)-t'^o{f^- p)}- 

0<a<s 0<p<M 

In the last step we substituted in (|2.5p . Formula (j2.30p follows from this by some calculus. 

From the definition (|2.22p of follow inequalities analogous to (|4.17p . and then 

with arguments like those in the proof of Lemma 14.51 we derive a right tail LDP 

(5.1) hm logP{zf2!L| > = -JeMri^) = - inf {RaDJs-a,t)(.r) 

n-i-oo L™*J:L™J ae[0,s] 

where Ra is the rate function from (|4.18p . For > the l.m.g.f. in (j2.24p satisfies 
~ '^ehor(^)- This would be a consequence of Varadhan's theorem if we had a 
full LDP, but now we have to justify this separately and we do so in Lemma |5. II below. 1 
Proceeding as in (j4.4ip and using ()4.50p 

A^r(;t)(o= sup iRm + j:.,M 

aG[0,s] 

= sup inf {aMe{0 + {s-a)Mp{C)-tM^,^^{-()}. 

ae[0,s] pG(C,At) 

Formula (|2.3ip follows from some calculus. The sup and inf can be interchanged by a 
minimax theorem (see for example |19) ) and this makes the calculus easier. □ 

Lemma 5.1. Let the partition function given by ()2.22p and let Je^hovi''') o,s given 

by dSl]). Then for Q < ^ < 9 

hm n-MogEe^'°^^Lnsj,Lntj = g^pjr^ - Je,hor(r)} = JehorlO- 

n->oo 

Proof. Let < ^ < 6*. Set 

7 = lim n log Ee L"sJ , and 7 = lim n log Ee L^sJ , L"tJ . 



n— >oo 



First an exponential Chebyshev argument for a lower bound: 

logP{log 2'f2'|'" > nr} < -^r + n'^ log Ee^ '"'^ ^L-Untj . 

Letting n — )■ 00 along a suitable subsequence gives ^ > — Jq hor('") for all r € M. Thus 
7>^e,hor(^)holds. 
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For the upper bound we claim that 

(5.2) hm lEI n-MogE(e^'°'^^L-tr"tJl|logZ[^^'''°^i > nrj) = -oo. 

Assume for a moment that (j5.2p holds. To establish the upper bound let (5 > and 
partition M with rj = i5, i (z Z. 

n-ilogE(e^'°'^^Lniuntj) <n-Mog[ e^^'^'+iPOog > nrJ 

(5.3) + e^'-- + E(e«'°s^Ln'u-J l{log zf^jl^Lj > ^^-^l 
By (|5.2p . for each M > there exists m = m(M) so that 

n-ilogE(e^'°'^^Lnlu-Ji{logzf2j% > nr„}) < -M. 
A limit along a suitable subsequence in (j5.3p yields 

7 < max {^ri+i - J0,hor(ri)} V ^r_m V (-M) 

— m<i<m 

< (supl^r - Je,hor{r)} + e<5) V er-„^ V (-M). 

The proof of the lemma follows by letting 6 —?■ 0, m — > oo and A4 —?■ oo. 
Now to show ()5.2p . Note that there exists a > 1 such that a.^ < 9, 



, „(9),hor \ 1/n 

Ee"«i°*5^L-J.Ln*J ) <oo. 

To see this, distinguish cases where < 1 or otherwise. Let N denote the number of paths 
and recall that N < e"^" for some c > 0: For < 1, 



For > 1, Jensen's inequality gives 

7(6), hor \ 1/n 



[nsj + [nt J 



/ ^{e),hov \ 1/n / r/ i-J: ^ \ afi \ 1/n / , i-^ ,N 1/" 

(Ee"^'°s^Ln.j,Ln*jy =(^E^(^ J2 n ]) -(^ n 



^6n(|^„3j_L„tj) 4=1 i=l 



< e""«M0(aC)*+^ 
To show (j5.2p . use Holder's inequality, 

n-HogE(^e^'°^^l^^TUl{logZl^;'^^^^^^ > nr} 

< a-Mogsup(Ee"^^°^^Lnluntj)i/« + _ l)a-^n-i logP{log zf^tLi > nr}. 

n L J 1 L J 

Taking a limit n — )• oo we conclude 

(5.5) logE(e^'°^^LniM-J l{log Z^^;"^"^^^^ > nr}) < Ci - C2Je,hor(r), 
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for positive constants Ci, C2. Letting r — >■ 00 finishes the proof because hnir^-oo >^0,hor('^) = 
00. ' □ 

Proof of Theorem \2.11[ We can assume 0<S,<9A{fi — 6) because otherwise the boundary 
variables alone force the l.m.g.f. to blow up. 

Let us record the counterpart of (I2.3ip for • Condition (j2.28p becomes 

(5.6) t(^o(M -0)- ^oif^ -0-O)> s{^o{0 + 0- ^o(^)). 

The conclusion becomes that the limit in (|2.25p exists and is given by 

.ver if dSSI) holds 

^ ^ ^ ''^'''^^^^ ~ \AtAO = AsAO, if m fails. 

The logarithmic limits lead to the formula 
(5-8) A,,(,,,)(0 = A^-,)(OvA^f^,)(0 

and we need to justify that this is the same as the maximum in (|2.20p . This comes from 
several observations. 

(i) As.t(^) = Jst{C) is always bounded above by the first branches of both (j2.3ip or ()5.7p . 
This is evident from equations ()4.49p - (j4.50p . 

(ii) Conditions (j2.28p and (|5.6p together define three ranges for (s,t): 

(a) (I2:28]l and 1^ both hold iff ait < s < a2t; 

(b) (I2:28]l holds and 1^ fails iff s > 02^; 

(c) (12:28]) fails and 1^ holds iff s < aii. 

The constants < ai < 02 can be read off ()2.28p and (j5.6p and the strict inequalities 
justified by the strict concavity of ^'o- 

(iii) In the maximum in (j2.20p we have 

(5.9) sMe{^) - tM^_e{-0 > tM^MO " sMe{-0 

iff s > a^t for a constant 03 > that can be read off from above. Strict concavity of ^'o 
implies that < ai < as < 02- 
Now we argue that 

(5.10) Ae,(s,t)iO = max{sMe(e) - tM^_e(-0, tM^-e{0 - sMe{-0}. 

This is clear in case (a) as this maximum is exactly Ag°^^^(^) V Aq^^^^^{^). In case (b), 
Ag™ (^) equals the left-hand side of (j5.9p which dominates both the right-hand side of 
(|5.9p and A^ ^(^). Consequently in case (b) also (|5.8p is the same as (|5.1Up . Case (c) is 
symmetric to (b). This completes the proof of (j5.10p . 

With one additional observation we can verify Remark I2.15[ Namely, As^tiO is in fact 
strictly bounded above by the first branch of either ()2.3ip or (j5.7p . The claim is easily 
verifiable when either of conditions (b), (c) are in effect. To see the strict domination 
when (a) holds, note that the unique minimizers in formulas (|4.49p ~ (j4.50p are linked by 
p = fj, + — 9. But if these formulas matched both first branches in (j2.3ip and (j5.7p . the 
connection would have to he p = fi — 9. This together with (|5.10p implies that As^((,^) < 
Ae,(.,i)(e) for all^G (0,;u). ' □ 
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